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Two parametrizations for second order velocity moments, the Batchelor 
O ■ parametrization for the r-space structure function and a common parametrization 

Q . for the energy spectrum, E(p) oc p~^/^exp(— p/p^^), are examined and compared. 

■^ ! In particular, we investigate corrections to the local scaling exponents induced by 

finite size effects. The behavior of local r- and p-space exponents differs dramati- 
cally. The Batchelor type parametrization leads to energy pileups in p-space at the 



^ • ends of the ISR. These bottleneck effects result in an extended r-space scaling range, 

Q . comparable to experimental ones for the same Taylor- Reynolds number Rex. Shear 



c^ ! effects are discussed in terms of (global) apparent scaling correction S(^'^pp to classical 

scaling. The scaling properties of 5C,°''^^{Re\) differ not only among the parametriza- 
tions considered, but also among r- and p-space for a given parametrization. The 
difference can be traced back to the subtleties of the crossovers in the velocity mo- 
ments. Our observations emphasize the need for more experimental information on 
crossovers between different subranges. 
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1 Introduction 

In the theory of fully developed turbulence, scaling ranges of velocity moments in r- 
and in p-space are often put into a one-to-one correspondence with each other. The 
two scaling exponents associated with the scaling ranges are believed to be equiva- 
lent. While this view is correct for an infinite system (with an infinite scaling range) 
the relation between r-space and p-space exponents becomes more complicated for 
finite Reynolds numbers. It is the aim of this paper to quantitatively examine these 
finite size effects, as they might well be essential to interpret experiments |jl|, ||, |^ 
and numerical simulations 0, |^ . 

We begin our short review of previous work relevant to the present investigation 
by defining the infinite scaling range exponents. In r-space, they are defined via 
the velocity structure functions 

DM(r) = {{u{x + r) - uix))"^) oc /-". (1.1) 

From a theoretical point of view [^, 0, Q , the p-space scaling exponents correspond- 
ing to the (discrete) Fourier transformation u{p) of u{x) are more easily accessible, 

{\u{p)r)ocp-'^'-. (1.2) 

Kolmogorov's classical dimensional analysis of the turbulence problem [0 gives of 
course the same result for both kinds of scaling exponents, namely (m = m/3. Since 
Landau's famous footnote in ref. |T^ it has been a matter of interest whether there 
are scaling corrections 6(rn = Cm — fri/S to the classical result in the limit of infinite 
Reynolds number i?e [|n], |T2], 0, g. 

Grossmann and Lohse discovered and analyzed finite size scaling corrections in 
their reduced wave vector set approximations (REWA, see [|r^, |T^| and refs. therein) 



of the Navier-Stokes equations. These calculations fill the huge gap in the Taylor- 
Reynolds number Re\ between experiments with Re\ > 10^ |TB, R[ and full numerical 



simulations, which reach only Re\ ~ 200 H . Similar to the highest Re\ experiments. 



REWA [|T4|, |T6| achieves Rex > 10'^, thus resolving four decades in p-space. Three 
different ranges were distinguished in p-space: The stirring subrange (SSR, small 
p) with slight intermittency (i.e., scaling corrections), the viscous subrange (VSR, 
large p) with strong intermittency, and the inertial subrange (ISR, medium p) with 
hardly any intermittency [|13|, |l^. The physical origin of the VSR and SSR scal- 



ing corrections was extensively discussed in ||T3. The VSR corrections arise from 



the competition between turbulent energy transfer downscale and viscous damping 



T7| , p!q] . The SSR scaling corrections are presumably due to the broken symmetry 
of the Navier-Stokes dynamics because of the finite size of the system: For small p 
only downscale energy transfer is possible, i.e., the translational invariance and the 
self similarity of the turbulent flow is broken by the boundary conditions. In addi- 
tion to the investigation of these local scaling corrections, REWA also offered the 
opportunity to study the i?eA-dependence of global corrections to classical scaling. 



Values of Re\ ranging from 10 to 10 could be simulated |T^, |19| and it was shown 
| TP|] that 5C,ra oc Re\~^'^ due to the spectral corrections to classical scaling, 

{\u{p)n^p-'^''\l + a^ij\ j, (1.3) 

which result from large scale anisotropy (e.g., shear). The p"^/^ shear correction 
has first been suggested by Lumley |^n|, who employed dimensional analysis, and 
was later also found in refs. pT] , p2| , |19|, ^ with the help of dimensional analysis in 



terms of Clebsch variables. In refs. [21, 241 the correction has been associated with 



conserved helicity flux in p-space. The parameter ps is the typical scale set by the 



strength of the shear s, ps = v-s^/e, and am a dimensionless parameter, presumably 
on the order of 1. It is not yet clear, whether the second term in eq.( |1.3|) is only a 
small-p correction or whether pure shear energy spectra E{p) oc (|tt(p)p)/p oc p~''l'^ 
exist, i.e., whether ps 3> Pl can be achieved. Here, pi = l/L is the momentum scale 
set by the external stirring force. If pure shear spectra exist, they will probably 



be more pronounced in cross spectra |25] £'12 (p). Yakhot [^] recently discussed 



experimental indications for pure shear spectra. Whether Maloy and Goldburg's 
p6[| scaling of the velocity structure function D^'^\r) oc r^/^ in Taylor-Couette flow 
with oscillatory inner cylinder corresponds to ( |1.3D is also not clear. 

A systematic analysis of the properties of r- and p-space scaling exponents in 
flnite--ReA turbulence has been performed by the present authors in ref.|27]. We 



demonstrated that the r-space crossover from the ISR to the VSR and from the 
ISR to a large-r saturation range can lead to energy-pileups at both ends of the 
p-space ISR, the so-called bottleneck phenomenon |^ . In other words, monotonous 
local r-space scaling exponents may give rise to non-monotonous local p-space expo- 
nents. Both physical (based on the conserved energy current in p-space) and formal 
explanations for this effect as well as a comparison with available experimental and 
numerical data were given in ref . [^ . 

In the present paper, we continue and extend our investigation of flnite size 
effects on local scaling exponents in r- and in p-space. In the absence of ana- 
lytical techniques, that would enable us to study these questions directly from a 
dynamical point of view (which would of course be preferable) , our strategy is the 
following: Throughout the paper we will compare two common parametrizations 
(and flnite size corrections thereof) for the scaling behavior of velocity moments, 
namely the p-space parametrization ( p.l| ) discussed by Foias, Manely and Sirovich 
p9|| (henceforth called 'FMS-Parametrization' for simplicity) and the Batchelor- 
parametrization (|2.2|) PO, |2H], common in r-space. The FMS-parametrization has 



already been discussed in earlier work by Brachet et al. |31] and Frisch et al. |18 
We will see that the main difference between FMS- and Batchelor-parametrization 
lies in the description of the crossover from the viscous to the inertial range. In 
this paper, these parametrizations are considered to be two examples for possible 
crossover scenarios and we hope that further experimental results will allow for an 



unambiguous discrimination between them. We believe that pointing out the impor- 
tance and the surprising consequences of the nature of the crossovers is a main novel 
aspect of our work, since up to now research was mainly concerned with infinite Re 
scaling exponents. 

The paper consists of two major parts: The first one deals with local small-p 



scaling corrections. The bottleneck energy pileup at the infrared end of the ISR ||2^ 



as well as the small p scaling corrections found in REWA |T3| and also by dimensional 



analysis in terms of Clebsch variables |jT9|, ^ ^] belong to this category. We 



believe that these are all manifestations of the broken Navier-Stokes symmetry due 



to the boundaries, i.e., large scale anisotropy [|T3l. To investigate possible relations 
among these effects we modify the above-mentioned parametrizations to describe 
the scaling corrections from REWA, examine the ensuing consequences in r-space, 
and perform a quantitative comparison for the local scaling exponents C2(p) resulting 



from the three different approaches in refs.|13, |T^, 2? . 

In the second part of the paper we focus our attention on the Taylor-Reynolds 
number dependences of the apparent (global) scaling corrections 8(^2'^ due to 
shear effects. While it is probably not too surprising that Batchelor- and FMS- 
parametrization lead to different i?eA-dependences of the apparent scaling correc- 
tions 6(2^^, our result, that both parametrizations give rise to different behaviour 
of 6(2^^ in momentum and coordinate space, respectively, is certainly unexpected. 
E.g., we find that the p-space result 6(2'^^''^ oc Rex~^ of ref.[jl9|, |^ corresponds to 
§(^^PP'^ oc Re\~^''^ in r-space. This observation can be viewed as yet more evidence 
for the fact that finite size effects can render r- and p-space exponents inequivalent. 

Specifically, our paper is organized as follows: In Sect. 2 we introduce the FMS- 
and Batchelor-parametrizations and review our earlier calculations [^ comparing 
both parametrizations in p- and in r-space. Sect. 3 contains the corresponding 
analysis for power spectra and correlation functions in the time domain [^. In 
Sect. 4 we investigate the consequences of the SSR scaling corrections found in REWA 
and compare the small-p scaling corrections of refs.[0, |I^, |2^. Sect. 5 is devoted to a 
detailed study of shear effects, in particular, to the apparent scaling corrections 6(2^^ 
induced in r- and in p-space. Finally, Sect. 6 contains our summary and conclusions. 



2 Batchelor— versus FMS— parametrization 

2.1 Definitions and Fourier transforms 



To describe the behavior of energy spectra E{p), the FMS-parametrization [^, |T^, 
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EpMs(p) = Eoe^/^p'"^/' exp {-p/p,) (2.1) 



has frequently been used to interpret experimental [^, |35 and numerical P, |T3|, |T 
data. Here, Eq is the p-space Kolmogorov constant, e the energy dissipation rate, 
and pd characterizes the crossover to the viscous range. On the other hand, measured 



structure functions D^'^\r) are well described by the Batchelor parametrization ^^, 
^,|T|,|3|,§B|,0g, 

erV(3z/) 



DS\r) 



^Hin^y"' 



(2.2) 



where u is the viscosity, r] = z/^/^/e^^^ the Kolniogorov length, and b = 
27r(4/3)ii^o/5 = 6.0 — 8.4 the experimentally determined |2^ r-space Kolmogorov 
constant. The generalisation of both parametrizations to (2 7^ 2/3 is straightfor- 
ward and was considered in ||2^. The essential aspects of our present work do not 
depend on the precise value of (2- Clearly, both parametrizations neglect the finite 
size of the system since they do not contain a scale for the external stirring force. 
Furthermore, the FMS-parametrization eq.( p.l|) does not contain any energy pilup 
(or bottleneck effect) [E3, |2^. The velocity structure function for a given energy 



spectrum can be calculated through the Fourier transformation [^ 

sin (pr) 



D^i 



AJ^ E{p){l- 



pr 



dp. 



(2.3) 



Inverting this equation, i.e., calculating the energy spectrum from a given structure 
function, requires a short discussion. Let us consider turbulence in a large but finite 
domain, so that E{p) ^ as p -^ 0. Then the term involving no trigonometric 
function on the rhs of eq.( p.3| ) is finite, A j^ E{p)dp = D^'^\oo) < cxd. Physically, 
this term corresponds to the total energy in the fluid. We can now straightforwardly 
invert eq.(|2.3D to give 






pr sin(j9r) 



D^'Hr)-D^^h 



00 



dr. 



(2.4) 



In the limit of infinite system size /^'•^•'(cxd) grows beyond all bounds, rendering 
eq.( |2.4| ) ill-defined at first sight. However, since J^ pr sm{pr)dr oc pS'{p), this 
affects only singular contributions at the origin which we may safely discard. We will 
therefore always use eq.( |2.4| ) with the understanding that D^'^\oo) = 0. Formally 
this means nothing more but that the Fourier transformation of a function will not 
change (apart from the ^-function) if the function is shifted by a constant. With the 
help of the transformation equations (|2.3| ) and (|2.4| ) we can calculate the structure 
function corresponding to the FMS-parametrization ( |2.1| ) and the energy spectrum 
associated with the Batchelor-parametrization 






r) 



4Eor(-2/3) f5 
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PdT 



5/6 



givmg 



(2.5) 



and 



Eji{p) 



— —P^d-j-^ 



1 
47r3z/ 



Pd 



d^ 
dp- 

-5/3^ 



exp{ipr'^x) 

(l + x2)2/3 



p'^'Ku/eiP) 



dx 



P 



"'K,/,{p) 



(2.6) 



respectively. Here, we introduced some abbreviations for simplicity: The dimension- 
less constant A has the value A = (9r(l/3))/(y2^22/35r(2/3)), p = pr'^, and K^ is 
the modified Bessel function of third kind I^H]. The crossover from the inertial to the 
viscous range is characterized by r^, p^ for the FMS-parametrization and by r^, p'^ 
for the Batchelor-parametrization. The large- and small-r limits of D^'^\r) are 
required to give D^'^\r) = 6(er)*^^/^) and D^'^\r) = er^/(3z/), respectively. Eq.( p.2|) 
is obviously designed to meet these constraints and comparing the asymptotic re- 
lations with the appropriate limits of eq.f ^TSf ) fixes Eq = 56(r(4/3))~^/27 = 1.74 
and 

p^ 1 = (106/27)^/^7/ ^ 2.34r/, (2.7) 

where we chose b = 8.4 [^. Now, r^ and r^ are defined by equating the asymptotic 
limits, er^/(3z^) = 6(er)^/^, so that we arrive at 

rd = ra = {3b)^/^ri^ 11.257]. (2.8) 

We note, however, that although r^ and r^ are the same (by definition), D^ (r) 

(2) (2^ (2) 

shows a sharper crossover from VSR to ISR than Z)pMs(r), i.e. D^ (r) > Dfms(^) 
for all r. This can be seen in Fig.l of ref.[|2^ where we compared D^ (r) with 
DpMsir). Finally, the p-space crossovers pd and p'^ are defined by the cutoff in the 
exponential decay of the spectrum for large p. Thus, pd is determined by eq.(pj.lD 
and, since K^{p = pr'^) oc p^^/'^exp(— pr^) for large argument, we have 

p',-' = r',^ 11.25r/. (2.9) 

Note that the naive expectation that (p-space crossover) x (r-space crossover) ^ 27r, 
holds in neither case. For eqs. ( |2.1| ) and ( p.5| ) we have raPd ~ 4.8 whereas for eqs. 
( p.2|) and ( |2.6| ) we have simply r'^p'^ = 1. 



2.2 Bottleneck phenomenon 

In [^] we showed that in contrast to eq.( |2.1| ) the parametrization eq.( |2.2| ) contains an 
important physical phenomenon, the bottleneck effect ||2^. This becomes apparent 



when comparing the energy spectra i?pMs(p) and E-q{p) in Fig.l. For small p both 
functions coincide. Around p ~ p^, however, an energy pileup in the crossover 
region of E-q{p) becomes noticable, leading to a non-monotonous logarithmic slope 
d\ogE-B/d\ogp. The local r- and p-space scaling exponents 

alogr dlogp 

of both the Batchelor- and the FMS-parametrization are plotted in the right part 
of Fig.2. The minimal local p-space scaling exponent of the spectrum Eb{p) is 
0.44, i.e., the scaling correction is an order of magnitude larger than the discussed 



intermittency corrections |2^, O]. This effect could explain, that the spectra in 



numerical simulations [^ |^, ^ are flatter than the classical expectation, rather 
than being steeper as one might expect from possible intermittency corrections. To 
quantify this effect, we have introduced an effective (global) p-space scaling exponent 
(^■^■^{Re\) in ref. [^ which only slowly approaches its infinite Rcx value. Very recent 
measurements of C,'^^^ {Re\) by Zocchi et al. [^ showed exactly the same behavior 
forC//(i?eA) m- 

Furthermore, in ref.|2^ we considered a straightforward generalization of the 
Batchelor-parametrization, 



Df{r) 



oc r 



'1 2 

r^ +r 



-2/3 



{L' + r'f'\ 



(2.11) 



which accounts for the crossover from the ISR to a large-r saturation range induced 
by the finite scale L = 1/pl set by the external stirring force. This second crossover 
might well not be universal, but our parametrization agrees reasonably well with 
available data [^ ^ ^. Calculating the corresponding spectrum we obtained for 
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(2.12) 



where p = p/pL and 1^2(0, 6, c,z) denotes a generalized hypergeometric function 
il]. The most prominent feature of this expression is a second bottleneck pileup 
at the infrared end of the p-space ISR. The local logarithmic slopes of ( p.ll|) and 
( p.l2|) are plotted in the left part of Fig. 2. Both bottlenecks have the same phys- 
ical origin, namely the broken symmetry due to finite size effects. The symmetry 
breaking scale is introduced by the stirring force and the finite size of the vessel, 
wind channel, or atmosphere |[l|, |^ at the infrared end of the spectrum, and by vis- 
cosity at the large-p end of the ISR. Formally, both bottleneck energy pileups orig- 
inate from the sharp r-space crossovers defined by the Batchelor-parametrization 
eq.(^) 



The physical explanation 



]T1|, ^ builds on the constant energy flux 
T{p) ~ pu{p) J dpidp2u{pi)u{p2)S{p + Pi +P2) downscale in p-space. For a detailed 
discussion of the bottleneck effect we refer to ref.[^. 



2.3 Higher order moments 

The preceding analysis cannot easily be extended to higher order structure functions. 
The connection between, say, the fourth order structure function D^'^^r) and the 
corresponding fourth moment of u{p), 

D^*\r) oc j 5{p^+P2 + p^ + p^){u{p^)u{p2)u{p^)u{p^)) 



(n(exp {ipj ■ r) - l))^^^^Jp^dp2dp^dp^, 



(2.13) 



is considerably more complicated than eq.( |2.3| ). Therefore we have to restrict our- 
selves to a few general remarks. 

Neglecting intermittency corrections we assume as a first approximation that 
Vr{x, t) = u{x + r,t) — u{x, t) and u{p, t) are Gaussian distributed, so that we may 
simply factorize higher moments (for even m), 

{\u{p)n ex (|^^b)H-/^ 
D^"'\r) oc (D^^\r))"'^\ (2.14) 



Note, however, that the above assumptions are not independent. In a completely 
homogeneous medium the second moments in p-space are local (i.e., {u*{p)u{p')) oc 
6{p — p')) and the second line in eq.( p.l4| ) is a direct consequence of the first one. 
Of course, the assumption of Gaussian factorization in eq. ( p.l4| ) is at variance with 
the fact that odd moments do not vanish (e.g., D^^\r) < for large r due to 
Kolmogorov's structure equation 



The similar looking factorizations in eq.( p.l4 ) lead to quite different results con- 



cerning the m-dependence of the crossovers rj'™ and pj^ (or pj™^ ) between VSR 
and ISR. These lengths are defined as above by matching the asymptotic behavior 
for large and for small r and by the cutoff in the exponentials, respectively. In 
r-space we get 

'"i ~ ''^d ''^'d ~ constant (2.15) 

for all m, which is in agreement with recent measurements |3^, while in p-space 

p^™) = 2pJV^ = 2pd/m (2.16) 

becomes smaller with increasing m. (The same relation holds for p'^). Eq. ( |2.16| ) 
has been numerically confirmed to a high precision [|l^ . Thus, for increasing m the 
ISR becomes smaller and smaller in p-space, whereas it remains invariant in r-space. 
Technically, this is due to the fact that we compare two power laws in r-space and 
a power law with an exponential in p-space. An intuitive understanding, we hope, 
is provided by the following remark: Raising D^'^\r) to some power smoothes the 
transition from VSR to ISR and consequently reduces the corresponding spectral 
strength at large values of p. 

3 Frequency spectra 

Instead of performing our analysis for the structure function D^'^\r) = D{r) (for 
simplicity, we drop the index 2) and the spectrum E{p) we can do the same for 
the longitudinal structure function Dj^'{r) = DL{r) and the longitudinal spectrum. 
Eiip) |2^, which are connected with the experimentally most easily accessible time 



structure function D{t) = {{ui{t + r) — ni(t))^) and its frequency power spectrum 



1 r°° 
P{u) = / dTcos{uJT)D{T) (3.1) 

71 Jo 



by Taylor's hypothesis [^ P5| . 

As Batchelor's parainetrization is also an excellent fit to the directly measured 
time structure function D{t) [^, |3^, ^, we want to give the corresponding fre- 



quency spectrum P{uj) for completeness. Again we restrict ourselves to classical 
scaling. 

First, for the VSR-ISR crossover 



t" 



"""' " (l + (r/r.)T/3 - P-^' 

we obtain from eq. ( p.l|) 

P{u) oc -u-'/'^K./.iu) + u'/^Ku/d^). (3.3) 

Here, r^ characterizes the VSR-ISR crossover and a) = c^r^j. The local logarithmic 
slopes ({t) and ({uj) of ( PT^ ) and ( ^.31 ) (defined as in ( [2.1(J[ )) are plotted in the 
right part of Fig.3. Similar to C{p), the local slope C(^) is again non-monotonous, 
reflecting the bottleneck phenomenon. But for P{uj) (and thus -E'i(p)) it is only half 
as large as above: The local exponent C('^) has a minimum of 0.56, compared to 



0.44 of C{p), found in ^j, see also Fig. 2. 

Next, for the ISR-large r saturation range we assume (again led by experimental 
data 11, 11,11), 

D{r)ocr'/'{l + ir/rLry'^\ (3.4) 

This second crossover is governed by the large eddy turnover time-scale tl- The 
corresponding spectrum reads 

r(5/6) / ^,5 3 3 ., , 20 ., c^.ll 5 5 ^,^\ 1., ^,4 3 .,, 



''^"^°^7HXWr^^6'2'2'") + 27^ ^^^(¥'2'2'"V-r^''^(3'''2'" 

(3.5) 
where u = ujtl/2. The local logarithmic slopes of ( p.4|) and ( p.5|) are plotted in the 
left part of Fig.3. The maximal value of C (a;) is 0.73 instead of 0.77 found for ({p) 
in 1^. Again, the bottleneck effect is smaller in cu-space than in p-space. 

Our Fig.3 has to be compared with the experimental local slope of P{uj), which 
we show in Fig. 4, taken from Praskovsky and Oncley's recent paper Pj. While 
the comparison is definitely not conclusive due to the experimental noise, there 
nevertheless seems to be a certain tendency towards the formation of energy pilups at 
both ends of the ISR scaling range. For a discussion of the quantitative discrepancy 
we refer to section 6. 



4 Small— 2? scaling corrections 

In this section we examine the consequences of the p-space SSR scaling corrections 
found by Grossmann and Lohse | jIB| , |n[ in their numerical analysis. To this end 

9 



we modify the energy spectra E-q{p) and Ep^s{p) calculated in section 2 to include 
these corrections at the infrared end of the ISR. Two major problems arise in the 
course of this procedure: First, we have to introduce the external stirring force scale 
into the p-space parametrizations to define a finite range for the scaling corrections. 
Second, we have to calculate the structure functions corresponding to the modified, 
more complicated energy spectra. In the following we explain how we deal with 
these problems, present our results, and compare our findings with other small-p 
corrections [0, |^ that are being discussed. 

4.1 Infrared cutoff 

Apart from our discussion of the small-p bottleneck we have considered only ideal 
turbulence in an infinite spatial domain up to now. Real turbulence is restricted 
to a finite range, i.e., there is a maximal length scale L or, equivalently, a small 
wave vector cutoff pi = 1/L. In order to discuss the small-p bottleneck in section 
2 the Batchelor parametrization was generalized to include this length scale. Here, 
we have to modify the parametrizations E-q{p) and Ej,^s{p) of the energy spectra 
accordingly. Therefore we multiply the spectra by {2/ti) a,Tctan{{p/pLY^^^). This 
amounts to imposing energy equipartition E{p) oc p"^ on the unforced wave vector 
modes |]2S|, ^, |T^ with p '^ Pl- Note, however, that the small-p behavior of the 



spectrum in Navier-Stokes dynamics has not been firmly established up to now 



1^ and we adopt energy equipartition for small p only as one of several possible 
scenarios. We will come back to this point at the end of section 4.2. As an immediate 
consequence of the cutoff the corresponding structure functions saturate for large r 
at the constant value D^'^^oo) ^ D^'^^L). 

The finite maximal length scale L allows for the introduction of the Reynolds 
number Re or alternatively the Taylor-Reynolds number Rex = Xui^rms/^-, where 
A = ui^rms / {,diUi)rms is the Taylor length and z/ the viscosity. Let us express pd and 
p'^ in terms of L and Re\. We have e = Ceu\^^g/L with q ~ 1 known from grid 
turbulence experiments |H3. Here, we neglect the i?eA-dependence of q for small 



Re\ 1^. Note, that q can directly be connected with the Kolmogorov constant 



h 1^5], |2^, ^. On the other hand e = 15z/(9iMi)^^g [^. Using these relations we 

finally get 77 = lh^f'^c~^LRe^ or, with q ~ 1 and the relations ( |2.7| ) and (|2.9|) for 

Pd, P'd^ 

pt ^ l^LRel''\ (4.1) 

p',' ^ SQLRel''\ (4.2) 

This establishes the connection between the length scales r^ (or r^) and L, and the 
Taylor-Reynolds number Re\. 
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4.2 REWA scaling corrections and structure functions 

In the reduced wave vector set approximations of the Navier-Stokes equation, Gross- 
mann and Lohse |T^ found deviations S(m{p) from classical scaling when locally fit- 
ting the spectra. These deviations occur - as already discussed in the introduction 
- only for small p (SSR intermittency) and for large p (VSR intermittency) , whereas 
no scaling corrections were found in the p-space ISR |]13|. To investigate how the 
p-space scaling corrections act in r-space, we model a spectrum according to the 



numerical results in |T3[ and numerically Fourier transform it into r-space. We will 
focus attention on p-space SSR scaling corrections. A short calculation reveals that 
for any function E{p) with — C(p) ~ 1 = dlogE/dlogp a modification defined by 

E{p) = E{p) ^P'^P/ (4.3) 

p° 

leads to a local exponent of the type 

This is exactly what we want: Assuming that ({p) = 2/3 we get the modified 
exponent (^{p) = 2/3 + 6 for small p until the infrared cutoff sets in. Furthermore 
the scheme is very flexible. We can introduce positive and negative corrections 6 
and let them become effective for p <^ pb or p ':$> pb depending on the sign of [3. We 
mention that this procedure, when applied to the exact small-r result D{r) oc r^ in 
order to extend D{r) to the ISR, immediately leads to the Batchelor parametrization 
(for /3 = -2and6 = 4/3). 

For our present purposes we choose SSR corrections between 6 = 0.0 — 0.04, in 
the range of intermittency corrections discussed in the literature |2^, |TI|, |I|. For 
large p we have ({p) = 2/3, (and for very large p the exponential cutoff is supposed 
to become effective). The crossover is determined by pb, we choose pb = lOpi — l5pL, 
i.e., allowing for about one decade of large scale intermittency corrections in p-space 
as found in [^. The parameter j3 determines the smoothness of the transition, we 
choose P = 2. The two spectra 

(2 I „2\<5/2 ^ / 2 I „2U/2 

EpMs(p) = E^^sip r'^^/ , EM = EM ^P' s (4-5) 

pO pO 

were each multiplied by the small p cutoff (2/7r) arctan((p/pj^)^^/^+'^) and then 
Fourier transformed. The Fourier transformations of E^ip) and £'pMs(p) (includ- 
ing the arctan-cutoff ) were performed numerically, employing a routine designed to 
cope with the strongly oscillating integrand. As the strong oscillations in eq.( |2.3| ) 
are exponentially damped by the asymptotic behavior of E{p), there are no serious 
numerical difficulties. The local scaling exponents for p- and r-space (defined as in 
eq. (|2.10|) ) are shown in Fig. 5a and 5b, respectively, for a Taylor- Reynolds number 
of Rex = 3000, which is in the range of typical experiments |1|. 
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Let us first discuss the results for tlie FMS parametrization. Witliout any small 
p corrections ( [4.5D we liave about one decade of more or less constant ({p) ~ 2/3, 
Fig. 5a. Tlie corresponding structure function DpMs{r), however, does not scale, see 
Fig. 5b, where no scaling range can be observed from the local exponent C(^)- This 
demonstrates that the transformation from p- to r-space is not completely local: A 
reasonably well-defined scaling range in p-space is mapped to a r-space curve with 
only very poor scaling properties or, perhaps more appropriately, with no scaling 
range at all. In the REWA calculations of ref. [|l6l a very similar behavior of ({r) 



has been found, which is not surprising, as the spectra are well parametrized by the 



FMS parametrization [jT3|, |T^, |T^. The poor scaling properties of OpMsi^) P^t, in 
our opinion, a question mark behind the FMS-parametrization, since experimental 
data for structure functions for the same Rex exhibit much better scaling. 

Only by introducing the small p scaling corrections ([4.5|) to EpMs{p), i.e., by 
making the scaling properties worse in p-space (Fig. 5a), we get improved (but still 
poor) scaling of the r-space structure function (Fig. 5b). As examples, we chose 6 = 
0.02, Pb/pL = 10, which is about what was found in the REWA calculations [|13], |T3 



and 6 = 0.04, Pb/pL = 15. Our result completely agrees with our findings of section 
2: Non-monotonous p-space scaling exponents lead to nicer scaling properties in 
r-space, i.e., the small p scaling correction ( [4.5|) can be interpreted as artificially 
introduced bottleneck energy pileup on the infrared end of the p-space ISR. 

The same analysis is performed for the Batchelor parametrization. Now, due to 
the large energy pileup at the ultraviolet end of the p-space ISR, the small-r scaling 
properties of Dsir) are improved considerably. Yet the p-space arctan infrared 
energy cutoff still corrupts r-space scaling properties for large r. Again, this effect 
can be partly compensated by local infrared p-space corrections of type ( |4.5| ). With 
these corrections, Dsir) shows better scaling properties in r-space (Fig.5b) which 
are now comparable to the experimentally realized scaling ranges of about 1.5 — 2 
decades for that Taylor- Reynolds number Re\ = 3000 [0]. 

To summarize: The simple arctan orexp cutoffs of the p-space ISR scaling range 
lead to unrealistically short r-space scaling ranges. Only the energy pileups at both 
ends of the p-space ISR lead to a realistic scaling range of the structure function, if 
compared with experiment |^. Batchelor parametrizations of the crossovers (|2.11| , 
|2.12|) include these energy pileups and give realistic scaling ranges for given Rex. Our 
findings also explain, why the r-space scaling found in REWA is worse than that in 
the p-spectra [|^, as the latter is quite well described by the FMS parametrization 

Coming back to section 4.1, we briefly mention that it is only the type of the 
ISR to large-r saturation crossover in ( |2.11| ) which determines the very small p ^ 
behavior for Eb{p). For /3 = 2 we obtain Eb{p) cx: p, for /? = 4 it is Eb{p) oc p^. 
It is not clear, which behavior is the more realistic one. For the Euler equations 
the latter can be proven |J3| to be correct and simply reflects energy equipartition. 



-"^The bottleneck pileup found in [Q is quantitatively smaller than that following from eq. (2.6) 
so that the FMS-parametrization is still an appropriate fit. 

12 



m 



but for the Navier-Stokes dynamics the situation might well be different [^ 
particular, as the thermal energy is orders of magnitude smaller than the energy of 
the large scale eddies. 

4.3 Comparison of infrared scaling corrections 

As discussed in the introduction, the large scale anisotropy (boundary, shear) leads 
to eq. ( p. .31) by dimensional analysis. We wonder how this result compares with the 
local slope C2(p) resulting from ( p. 121) and from the REWA calculations [p!3| . 



Of course we can compare only the crossover region, as in the original work by 
Lumley [^ and also in the derivation in [^] the second term (oc p'"^^^) in eq.( |1.3D 
is considered to be small. 

The comparison for the local slopes C^ip) is given in Fig. 6. Clearly, although all 
curves show 5C,2{j>) = C2(p) — 2/3 > at the infrared end of the p-space ISR, they 
do not agree quantitatively. Note however, that no wave vectors smaller than the 
forcing scale were included in REWA. The presence of such wave vectors leads to the 



infrared bottleneck effect ||2^. It would be interesting to include such wave vectors 
in full numerical simulations or in REWA type simulations and to examine, whether 
a bottleneck energy pileup as in (|2.12|) or in (|3.5|) will occur. Also note from Fig.6, 



that the REWA corrections, which we introduced into the energy spectra by means 
of the parametrization chosen in the preceding subsection, could also be described 
by eq. (|1.3|) with 0^2 ~ 0.1. This confirms and justifies the analysis done in ref. 
|T^, see also section 5. The bottleneck corrections are comparable in size with the 



corrections due to eq. (|1.3|) (with the arctan cutoff) for a2 = 0.5. 

The quantitative discrepancies in Fig. 6 should not be too surprising. The small 
p spectral shape is far from being universal due to different boundary conditions 
and different kinds of stirring. Yet we believe that the three kinds of discussed 
small p scaling corrections 5(^2 (p) > all have the same origin, namely the broken 
symmetry of the Navier-Stokes dynamics for small p and the broken self similarity 
of the turbulent flow. 

4.4 Higher order moments 

In principle, the analysis of section 4.2. can be repeated for higher order velocity 
moments. In the REWA calculation [l^ it was found, that for large Rex higher 



order moments nearly factorize into second order moments, of course apart from 
the p-space SSR and VSR intermittency corrections. One could assume such a 
factorization and parametrize the small p scaling corrections by 

with Sm < for m > 3. We have performed this calculation but refrain from 
discussing the outcome in detail because nothing essentially new can be learnt. The 
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results simply confirm the views already developed when dealing with the second 
moments. 

5 Shear effects 

In this section we investigate the consequences of an extended shear scaling range. 
We will modify the Batchelor- and the FMS-parametrization in such a way that a 
crossover to a shear scaling range occurs at the scale of the stirring force pi = 1/L. 
We then proceed to calculate the impact of this modification on scaling exponents 
in both r- and p-space. In particular, we focus attention on the apparent scaling 
correction 5C,'^^^{Re\) [^ as a function of the Taylor-Reynolds number Re\. This 
correction is induced by the crossover itself and does not depend on the extension 
of the shear range. 

5.1 Shear parametrizations 

The generalized FMS-parametrization (cf. eqs. (|1.3|) and ( p.l|) ) 



{\u{p)n (X p-"^'' \i + a„A^\ 1 exp (-p/p,). (5.1) 




was shown |T^ to lead to an apparent scaling correction, defined by 

5Cr'" = min,(C^(p))-m/3, (5.2) 

where, as usual, Cmip) = —d\og{\u{p)\"^) /dlogp. For pi <^ Pd it was found that 

SCr^' = sign(a^)^ [^fc^U' = c^Re-'''' = c'^Ref^', (5.3) 

i.e., the apparent scaling corrections vanish with increasing Rex with a —3/5 power 
law for all m. Eq. ( [5.3| ) has been numerically confirmed by reduced wave vector 
set calculations |jT9|. The dimensionless constants Cm, c^, and am are found to be 
negative for m > 3 [|19| . This means that higher order moments do not factorize into 



second order moments in agreement with the p-space SSR scaling corrections |T^ 
found numerically. Note that if one assumes factorization of higher order moments 
in the shear range, shear can not account for experimentally observed [|l| scaling cor- 
rections. The prediction ( |5.3| ) ||TP[ is also in agreement with very recent experiments 



which clearly show a decrease of 6^^^'^ for increasing Rex [^ |. 

The local scaling exponent of eq. (|5.1|) is shown in Fig. 7a for Rex = 1500 (which 
is determined according to eq. ( |4.1j )). The apparent scaling correction dQ^'^ ~ 0.06 

is very large. Note, however, that according to eq.( p73D 6(2^^^'^ is proportional to 

3/5 
0:2 and will therfore be smaller for a2 < I- The Rex dependence is displayed in 

Fig. 8. For large Rex the asymptotic result ( |5.3| ) is recovered. 

14 



We have mentioned in the introduction that the form of the shear correction 
in eg. ( p. .31 ) is based on dimensional analysis only. One could argue that similar 



considerations as in ||20| can as well be directly applied to r-space expressions. 
Assuming the transition to the shear range (with r-space exponent 4/3) to be of 
the Batchelor type we are led to the ansatz 

£7 rp ' rp ^ j "J 

""^ ' ?>u (r^f^ ^2)2/3 (^2 + ^2) -1/3' W-^^ 

with Ts = p~^ = Je/s^. We assume again that shear sets in at the stirring scale, 
hence r^ = L. For simplicity, we again dropped the index 2 of the structure function, 
as we will restrict ourselves to D^'^\r) from here on. Eq. (|5.4|) should be viewed 
as a generalization of the Batchelor-parametrization in complete analogy to the 
generalized FMS-parametrization eq.( |5.1j ). The local scaling exponent of Db{t) is 

given by 

o »,2 4- r^ 

CW = 2+777^—^-77^^—^, (5.5) 

and is shown in Fig. 7b. The apparent scaling correction is defined analogously to 



eq-dl) 



^^app,r ^ inin,(C(r)) - m/3. (5.6) 

For the same Re\ = 1500 as above, dC,""^^'"^ ~ 0.0028 is now much smaller than the 
corresponding value for the generalized FMS parametrization. This reflects the much 
better scaling properties of the Batchelor parametrization (compared to FMS, cf. 
fig. 5), which we have extensively discussed in sections 2 and 4. The Re\ dependence 
of (J^'^PP''' is displayed in Fig. 8. For large L^ r'^ (i.e., for large Re\) we obtain 

^^app,r ^ ^2-^ oc i?e-3/4 qc i^eA^'/^ (5.7) 

3 ij 

Hence, also the asymptotic Rex dependence of 3C,"'^^''^ for the Batchelor parametriza- 
tion is quite different from that of (5(^"pP'P for the FMS parametrization. One might 
argue that this is not very astonishing since we compare two different parametriza- 
tions. Let us therefore Fourier transform both parametrizations and reexamine their 
scaling properties thereafter. 

5.2 Fourier transforms 

Transforming the FMS-parametrization eq.( |5.1| ) with the help of eq. ( PT^ ) to r-space 
leads to 

^ , , r(-2/3) /5 /^ 2 2\5/6 . /5 , , 

D^^s{r) oc -— - — ^ -Par - [1 + p^r \ sm - arctan(prfr) 
{5/3)rpJ Vo ^ ^ \6 

a2T(-A/3)p'^/^ /7 r 2 2\V6 . /7 , x i ^ /. on 

+ 7/3 [i^Vd-r - (1 + Pd^ ) sm ( - arctan(prfr) ) ) . (5.8) 
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The local slope of eg. ( ^.81 ) for Re\ = 1500, 0^2 = 1, is also plotted in Fig. 7a, 
to compare it with the slope of eq. ( |5.1| ). Now (5^"^^'^ ~ 0.11 is even larger than 
§(^cipp,p (^both for FMS), which is clearly understandable from Fig. 5 because of the 
even worse scaling properties of the FMS parametrization in r-space (compared to 
p-space) . 

What is more surprising is, that now for the same (FMS) parametrization the 
Re\ dependence of S^""^^'"^, cf. Fig. 8, is different from that of 5(^"pP'P. For the former 
we obtain the asymptotic result 

5C"PP'^(i?eA) oc Rex-'/\ (5.9) 



which is a considerable flatter dependence than (|5 

Finally, we calculate the spectrum corresponding to the generalized Batchelor 
parametrization ( ^.4|) . For r ^ r^ we can derive an analytical result, which we give 
in appendix A for completeness. In the general case, we have to restrict ourselves to 
a numerical treatment. The numerical Fourier transformation (|2.4|) of ( ^.4|) (or of 
( p.2|) ) is more delicate than the inverse transformation (|2.3|) due to the absence of an 
exponential cutoff. In appendix B we explain how we convert the strongly oscillating 
integral to a rapidly converging one by means of contour integration techniques |^ 



The result for the local p-space slope of the generalized Batchelor parametrization 
is shown in Fig. 7b. Of course it shows the ultraviolet bottleneck energy 
pileup 1^, which we had discussed in detail in ref. |^. But now, in addition. 



the spectrum shows reduced spectral strength at the infrared end of the p-space 
ISR, i.e., a decreased local slope ({p) < 2/3. This effect can be interpreted as, 
so to say, an inverse bottleneck effect and can both formally and physically be 
interpreted along the same line of arguments as the bottleneck pileups discussed 
above and in [^. Formally it reflects the sharp crossover from r^/^ to r^/^ scaling 
in the structure function. Physically fl^, ^, the constant energy flux T(p) ~ 
pu{p) J dpidp2u{pi)u{p2)S{p + Pi + P2) downscale now requires reduced spectral 
strength at the infrared end of the ISR, as the spectral strength is increased in the 
shear range. So it is just the opposite situation as that discussed in section 2.2, see 



also Ihs of Fig. 1 and ref. [^ . Correspondingly, there is also an energy pileup at the 



high p end of the shear range, which may be a consequence of the constant helicity 



flux in this region [^. It leads to a local slope ({p) > 4/3. 

It is not our primary goal to speculate about the nature of this crossover. What 
is important here for the discussion of apparent scaling corrections is, that positive 
local scaling corrections S({r) > in r-space lead to negative S({p) < in the p- 
space ISR, see Fig. 7b. Thus the apparent scaling corrections SC^^^'^ have to be 
defined as 

S(^PP'P = maXp^«p«p^(C2(p)) - m/3, (5.10) 

The i?e;^-dependence of 5^"^^'^ is shown in Fig. 8. From the data we conclude that 

SC^P'P oc -Rex'^^^ < 0. (5.11) 
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For Re\ = 1500 we have SC^^P'^ = —0.0082. This means that the quahtative differ- 
ence between the apparent scahng corrections in r- and in p-space is even larger for 
the Batchelor- than aheady for the FMS-parametrization. For the former, even the 
sign of the corrections is reversed as we go from r- to p-space, while for the latter 
only the magnitude and the asymptotic Rex scaling exponent change. 

Let us summarize the analysis of this section. We found that the subtleties of the 
ISR to VSR and ISR to shear range crossovers govern the Rex dependence of 6("'^^. 
Moreover, for both parametrizations discussed (generalized FMS and Batchelor) the 
scaling corrections are quite different in r- and in p-space. So, when analyzing the 
experimental /2eA-dependence of S^""^^, one should also expect different results in r- 
and in p-space, or, correspondingly, in the r- and u- domain. 

Note, that because of the considerable lack of experimental information about the 
shear range crossover, we consider eqs. (|5.1|) and (|5.4|) only as examples. As pointed 
out above, dimensional analysis cannot distinguish between them. We think it is 
worth while to experimentally or numerically study the crossover between ISR and 
shear range and hope that our analyses stimulate to do so. This might lead to a 
better understanding of the Rex dependence of scaling corrections. 

6 Summary and conclusions 

Throughout the paper we have demonstrated, that scaling properties in r- and in 
p-space can be quite different. One could argue that in the infinite Rex limit these 
differences are irrelevant. This is of course correct. Yet, as we demonstrated, for 
those Rex which can be achieved in experiments and even more so for the numerical 
ones, the finite size corrections are considerable and it is important to know what 
their influence is to be able to interpret the data correctly. Moreover, the apparent 
scaling correction S("-pp due to shear corrections even show asymptotically different 
Rex scaling behavior. For the FMS type parametrization we had obtained SC^^^'^ oc 
Rex~^^^ ^ 0] and ^C"^^'"' oc Rex~^^^, for Batchelor type parametrizations f^C^^'"' oc 



Rex~^^^ and -SC'pP'P oc Rex'^''^. 

Comparison of the size of the scaling ranges for given Rex between experiment [|I| 
and our parametrizations lets us favor a Batchelor type parametrization rather than 
a parametrization of FMS type. The latter, consisting of a power law in p-space 
with a large-p exponential cutoff and a small-p arctan cutoff, does not exhibit any 
bottleneck energy pileups at the ends of the p-space ISR and leads to unrealistic 
short scaling ranges in the r-space structure function. In other words, combining 
all regimes discussed in this paper (VSR, ISR, shear range and large-r saturation 
range) we think that the p-space parametrization 

E{p) = ?^2£!^arctan f ( ^] ] p-'/' ( I + ( ^] \ exp {-p/p,), (6.1) 
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with Pd^ Ps ^ Pl is less favorable than an r-space parametrization 

^ /4/3 -2/3 r-4/3 

^ ' 'iv {r'j + r2)2/3 (r2 + r2)-i/3 (L2 + r2)-2/3 ' ^''- ^ 

with r'^ < Tg < L, which shows bottleneck effects in p-space. In many isotropic 
turbulence experiments Ts ^ L and the shear range will be suppressed. If less 
isotropy is achieved in experiments, we may have, say, r^ k. L/4. In this case, the 
energy pileup due to large-r saturation (section 2 and ref. ||2^) and the spectral 
strength reduction due to shear effects will partly compensate each other at the 
infrared end of the p-space ISR. This leads to a smaller change of the local slope 
than predicted by eq. (|2.12|) or (|3.5| ). 

Another point to be kept in mind when comparing our predictions ( |2.12| ) and 
( p.5|) with experimental data is the issue of averaging. While experimental scaling 
exponents necessarily represent data averaged over a certain interval, we defined a 
pointwise (local) slope in ( |2.1CI| ) and used this basic quantity throughout the paper. 



In the numerical REWA calculations [|T^, |T4[, locally averaged instead of pointwise 



slopes were determined by fitting the parametrization ( p.l| ) to the data in each in- 
terval [p/a/IO, pv^IO]. To estimate the effect of averaging, we performed a running 
average of the pointwise slopes in fig. 3 (the quantity most easily accessible in ex- 
periments) using an averaging range of [a;/\/lO, a;\/TO], see fig. 9. As expected, the 
bottleneck pileups become attenuated, so that they are now quantitatively closer to 
the measured ones in fig. 4. Similar averaging algorithms could also be applied to 
the other local slopes. 

Finally, we mention that while eqs. (|6.1|) and ( |6.2|) explicitly distinguish between 
shear effects and large r saturation effects, caused by the boundary conditions (i.e., 
large scale anisotropy), these effects might be more interwoven. In section 4.3 we 
had compared shear effects, finite size effects, and the numerical REWA results. 
Qualitatively, they all lead to S({p) > for small p, but the quantitative agreement 
was less satisfactory. We again stress the necessity to produce as clean shear ranges 
as possible in experiments. Only then, one will be able to experimentally study the 
crossover phenomena associated with the shear range, which - together with those 
from VSR to ISR - might well be a key in understanding scaling corrections. 
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A Analytical Batchelor— type shear spectrum 



We can analytically perform the transformation of the Batchelor-parametrization 
eq. (|5.4| ) under the assumption that r 3> r^, so that 






(A.i; 



With this approximation we get 



Eip) 



1 r°° 
~27r Jo 
p e 



prD^'^'{r) sin(pr)(ir 



271 'iurl'^ dp \ 
The integral in eq.( |A.2| ) can be solved |^ to give 



r'^'^ir'l + r^) ^'^ cos{pr)dr. 



(A.2) 



'^ ('"s + '"^) ^ cos{pr)dr 



' 12/3 ^6' 2' 6' 4 ^ 2 ^ ' '^ ^3' 3' 6' 4 ^' ^ ^ 



where p = pvs- Reinserting this result into eq.( [A.2| ) we derive after some algebra 
the following, rather clumsy expression, 



^, , erf/3 I 15 r(5/g)2 
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For small p the local scaling exponent of E{p) is 4/3, for large p it is 2/3, but the 
transition from one range to the other is non monotonous. This is reflected in the 
right part of Fig. 7b. In the left part of that figure, the ultraviolet bottleneck energy 
pileup can be seen in addition, which is not included in eqs. ( |A.1| ) and ( [A.4| ). 



B Contour integration for oscillating integrands 

The numerical Fourier transformation ( p^.4| ) of ( p.4| ) cannot straightforwardly be 
performed, as the integrand is strongly oscillating and not exponentially damped. 
To cope with this problem, we employ contour integration techniques ESI. Plugging 



4]) into (|2.4| ) we obtain after some algebra 



^W- 3^^d r, 4^^^3y_^ (^.2 + ^2)2/3 ^^- ^^■') 

The integral has singularities or zeros at ±ir^ and ±irs. Taking the correct branch 
cuts and performing the corresponding contour integration in the upper half plane, 
we obtain 

oo (r^ + r^)^/^exp (ipr) ^ r°° \z'^ — r'^\^/^ exp {—pz) m o\ 

-oo {r2 + r2)2/3 J^,^ (^2 _ r^2-)2/3 

or 

E{p) = — r/ r, '■'-— — ^-— dz, B.3 

Su Att Jr'^ [z^ — rJi) '^ 

which can now be straightforwardly integrated. Our numerical result is displayed in 
Fig.7b. 
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Figure Captions 



Figure 1: In the right part of the main curve we show the energy spectrum eq. ( p.6|) 
(solid) with, and the spectrum eq. ( pTTP (dashed) without the energy pileup. In the 
left part the spectrum due to ( |2.12|) is shown. In the insets, the spectrum is enlarged 
around the energy pileups and compared to classical — 5/3-scaling. 



Figure 2: The local p-space scaling exponents ({p) (solid), and the local r-space 
scaling exponent C('" = 1/p)- 



Figure 3: The local r- and cu-space deviation from classical scaling, S({t = 1/uj) 
and S({uj), when assuming Batchelor kind crossovers (|3.2| ) and 



Figure 4: The experimental local p-space deviation from classical scaling, S({p), for 
the longitudinal energy spectrum Ei{p). This curve corresponds to the experimental 
S({u!) via the Taylor hypothesis. The data are taken from Praskovsky and Oncley 
H with kind permission of the authors. 



Figure 5: (a) Local p-space scaling exponents ({p) of the FMS and the Batchelor 
type energy spectra, both with the arctan cutoff for small p, see text. To allow for 
comparison with r-space. Fig. 5b, we plotted ({p) versus Pl/p rather than versus 
p/Pl- The Taylor Reynolds number is Rex = 3000, cf. eqs. (|4.1j , [4.2|) . From bottom 
to top on the rhs of the figure, the three pairs of curves correspond to (i) no small p 
scaling corrections, (ii) small p scaling corrections according to ([4.5|) with 6 = 0.02 
and Pb/pL = 10, and (iii) small p scaling corrections with 6 = 0.04, Pb/pL = 15. 
(b) Local r-space scaling exponents ({r) for the six curves of Fig. 5a. From bottom to 
top for both the FMS and Batchelor triple of curves: no scaling correcions, 6 = 0.02, 
Pb/PL = 10 and 6 = 0.04, ph/pi = 15. 



Figure 6: Local scahng corrections 6(2{,p) = C2{p) — 2/3 due to the REWA calcu- 
lations [l^, Q, due to the infrared bottleneck formula ( p.l2| ), and due to eq. ( |1.3|) 



(with the arctan cutoff for small p) with three different values for the unknown 
parameter a2- 
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Figure 7: (a) Local scaling exponents ({r) (solid) and ({p = 7/r) (dashed) when 
shear corrections according to ( p.l\ ) are present, with Rex = 1500 {pd = 3227) and 
^2 = 1- The parameter 7 serves to shift the p-space curve slightly to ensure that 
the minima of the two curves coincide. 

(b) As in (a), but now shear corrections according to ( |5.4| ) and its Fourier transform. 
Again, we chose Rex = 1500, {p'^ = r^^ = 676). 



Figure 8: Double logarithmic plot oi dC^^'^iRex) and 6CPP'P{Rex) for FMS (as = 1) 
and Batchelor parametrization. 



Figure 9: Same as in fig. 3, but now in addition the averaged local slopes. The 
averaging range is [iu/^/TO,uj\/TO]. 
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* On leave of absence from Fachbereich Physik, Universitat Marburg, Renthof 6, 
D-35032 Marburg. 
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